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DIRICHLET-NEUMANN WAVEFORM RELAXATION METHOD 
FOR THE ID AND 2D HEAT AND WAVE EQUATIONS IN 
MULTIPLE SUBDOMAINS 

MARTIN J. GANDER*, FELIX KWOKt, AND BANKIM C. MANDAL * 


Abstract. We present a Waveform Relaxation (WR) version of the Dirichlet-Neumann algo- 
rithm, formulated specially for multiple subdomains splitting for general parabolic and hyperbolic 
problems. This method is based on a non-overlapping spatial domain decomposition, and the it¬ 
eration involves subdomain solves in space-time with corresponding interface condition, and finally 
organize an exchange of information between neighboring subdomains. Using a Fourier-Laplace 
transform argument, for a particular relaxation parameter, we present convergence analysis of the 
algorithm for the heat and wave equations. We prove superlinear convergence for finite time window 
in case of the heat equation, and finite step convergence for the wave equation. The convergence 
behavior however depends on the size of the subdomains and the time window length on which the 
algorithm is employed. We illustrate the performance of the algorithm with numerical results, and 
show a comparison with classical and optimized Schwarz WR methods. 


Key words. Dirichlet-Neumann, Waveform Relaxation, Heat equation, Wave equation, Domain 
decomposition. 


1. Introduction. A recent version of Waveform Relaxation (WR) methods, 
namely Dirichlet-Neumann Waveform Relaxation (DNWR) has been introduced in 
[16, 10, 9] to solve space-time problems in parallel computer. This iterative method 
is based on a non-overlapping domain decomoposition in space, and the iteration re¬ 
quires subdomain solves with Dirichlet boundary conditions followed by subdomain 
solves with Neumann boundary conditions. For a two-subdomain decomposition, we 
have proved superlinear convergence for ID heat equation, and finite step convergence 
for ID wave equation. In this paper, we extend the DNWR method to multiple subdo¬ 
mains, and present convergence analysis for one dimensional heat and wave equation. 
We also present convergence result for two dimensional wave equation. 

In a different viewpoint, the WR-type methods can be seen as an extension of 
DD methods for elliptic PDEs. The systematic extension of the classical Schwarz 
method to time-dependent parabolic problems was started in [11, 1 ]; later optimized 
SWR methods have been introduced to achieve faster convergence or convergence 
with no overlap, see [ ] for parabolic problems, and [ ] for hyperbolic problems. Re¬ 
cently Neumann-Neumann Waveform Relaxation (NNWR) algorithm is formulated 
from substructuring-type Neumann-Neumann algorithm [!, 26, . ] to solve space-time 
problems; for more details see [13, 14, 15], The DNWR method thus can be regarded 
as an extension of Dirichlet-Neumann (DN) method for solving elliptic problems. The 
DN algorithm was first considered by Bjprstad & Widlund [ ] and further studied in 
[3], [ ] and [ ]. The performance of the algorithm is now well understood for elliptic 
problems, see for example the book [21 ] and the references therein. 

We consider the following two PDEs on a bounded domain fl C R d ,0 < t < T, 
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d = 1, 2, 3, with a smooth boundary as our guiding examples: the parabolic problem 


du 

~dt 
u(x , 0) 

u{x , t) 


V ■ (k(x, t)Vit) + f(x,t), 

x G H, 0 < t < T, 

u 0 (x), 

iGfl, 


x G <9H, 0 < t <T, 


where n(x,t) > 0, and the hyperbolic problem 


c 2 (x)Au 

= 

x G H, 0 < t < T, 

u(x , 0) 

= v 0 {x), 

a)GH, 

u t (x, 0) 

= w 0 (x), 

a)GH, 

u(x, t) 

II 

x G (9H, 0 < t < T, 


( 1 . 1 ) 


( 1 . 2 ) 


with c[x ) being a positive function. 

We introduce in Section 2 the non-overlapping DNWR algorithm with multiple 
subdomains for (1.1), and then analyze its convergence for the one dimensional heat 
equation. In Section 3 we define this method to hyperbolic problems, and analyze 
convergence behavior for one dimensional wave equation. We extend our result to two 
dimensional wave equation and prove similar convergence behavior as in ID in Section 
6. Finally we present numerical results in Section 7, which illustrate our analysis. 


2. DNWR for parabolic problems. The Dirichlet-Neumann Waveform Re¬ 
laxation (DNWR) method for parabolic problems with two subdomains is introduced 
in [10, 16]. In this section we generalize the algorithm to multiple subdomains in one 
spatial dimension. We present different possible arrangements (in terms of placing 
Dirichlet and Neumann boundary conditions) and with a numerical implementation 
of these arrangements for a model problem we determine the best possible one. We 
then formally define the DNWR method for (1.1), and analyze its convergence for the 
one dimensional heat equation. 


2.1. Motivation. Suppose we want to solve the ID heat equation 


du 

~dt 
u(x, 0) 

u(x, t) 


A u, x G H, 0 < t < T, 

uo(x), x g n, 

g{x, t ), x G 5H, 0 < t < T, 


( 2 . 1 ) 


using the DNWR method. The spatial domain H = (0, 5) is decomposed into five 
non-overlapping subdomains H^ = (xi-±,Xi),i = 1,..., 5, see the left panel of Figure 
2.1, with three possible combinations of boundary conditions along the interfaces, 
right panel of Figure 2.1 and two arrangements in Figure 2.2. D in blue denotes the 
Dirichlet condition along the two physical boundaries, whereas D and N in red denote 
the Dirichlet and Neumann boundary conditions along the interfaces. We are given 
Dirichlet traces {g?(£)}. = as initial guesses along the interfaces {xi}f =1 . 

First arrangement (Al):. Here we extend the two subdomain-formulation [10, 
] to many subdomains in a natural way, see the right panel of Figure 2.1. With the 
intial guesses, a Dirichlet subproblem is solved in the first subdomain Hi, followed by 
a series of mixed Neumann-Dirichlet subproblem solves in the subsequent subdomains 
(H^ i = 2,..., 5), exactly like in the two-subdomain case. Thus the DNWR algorithm 
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Fig. 2.1: Decomposition of the domain on the left, and arrangement of boundary 
conditions due to Al on the right 



Fig. 2.2: Two different arrangements of boundary conditions: A2 on the left, and A3 
on the right 


is given by: for k = 1, 2... and for 9 £ (0,1] compute 


d t u\ - d xx u\ 
u\(x, 0) 
«i(0 ,t) 
Ui(x!,t) 

and for i = 2,..., 5 

dt'Ui - d xx u\ = 

Ui(x, 0) = 

d x Ui(xi-i,t) = 
u i {Xi , t) = 


= 0, in Di, t > 0, 

= u 0 (x), in Di, 

= g(o,t), t> o, 

= t> 0, 


0, 

Uo(x ), 

at 1 ®, 


in tti, t > 0, 
in f li, 
t > 0, 
t > 0, 


with (t) = g( 5, t) for the last subdomain along the physical boundary. The updated 
interface values for the next step are then defined as 

9i®=8ul 1 (x it t) + (l-0)gt 1 ®. 

We now discretize (2.1) using standard centered finite differences in space and back¬ 
ward Euler in time, and solve the equation numerically using the above algorithm for 
different time windows. For the test we choose iio(x) = 0,g(x,t) = (x + l)t,gf(t) = 
t 2 ,t £ [0,T]. Figure 2.3 gives the convergence curves for different values of the pa¬ 
rameter 9 for T = 2 on the left, and T = 20 on the right. 
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Fig. 2.3: Arrangement Al: Convergence of multi-subdomain DNWR for various re¬ 
laxation parameters 9 1 on the left for T = 2 and on the right for T = 20 


Second arrangement (A2):. This is the well-known red-black block formula¬ 
tion, described in the left panel of Figure 2.2. In this arrangement, we solve a Dirichlet 
subproblem and a Neumann subproblem in alternating fashion. Given initial Dirichlet 
traces along the interfaces, a series of Dirichlet subproblems is first solved in parallel 
in alternating subdomains (12i, SA 3 , ST 5 ), and then a series of Neumann subproblems 
is solved in the remaining subdomains (£^ 2 j ^ 4 )- So this type of DNWR algorithm is 
given by: for k = 1,2 ... compute for i = 1,3, 5 


dtu’l - d xx u\ 
u i( x > 0 ) 

Ui(Xi-l,t) 


0 , in f 2j,i > 0 , 


u 0 (x), in 12 ,, 

9i-i(t), t> 0, 
9i 1 {t ), t > 0, 


with 3 §(f) = g(0,t),g£(t) = g(5,t), and for * = 2,4 


d t u\ - d xx u% 
u-(x,0) 
d x Uj(xi-i,t) 
d x u^(xi,t) 


0 , 

Uo(x), 

d x Ui_i(xi-i,t), 

d x u^ +1 {xi,t), 


in fli, t > 0 , 
in 12 i , 
t > 0 , 
t > 0 , 


together with the updating conditions 

9i(t) = eu k i+l (xi,t) + (1 - 0)g k ~ 1 (t), i = 1,3, 
9i(t) = 0u k (xi,t) + (1 - i = 2,4, 


where 9 £ (0,1] is a relaxation parameter. 

We now implement this version of DNWR algorithm for different time windows, 
picking the same problem and initial guesses as for Al. Figure 2.4 gives the conver¬ 
gence curves for different values of the parameter 9 for T = 2 on the left, and T = 20 
on the right. 

Third arrangement (A3):. We now consider a completely different type of 
arrangement, proposed in [4] and shown in the right panel of Figure 2.2. Given initial 
guesses along the interfaces, we begin with a Dirichlet solve in the middle subdomain 
D 3 , followed by mixed Neumann-Dirichlet subproblem solves in the adjacent subdo¬ 
mains, in an order D 2,124 first and then in 12i, G 5 . This third version of the DNWR 
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Fig. 2.4: Arrangement A2: Convergence of multi-subdomain DNWR for various re¬ 
laxation parameters 9 1 on the left for T = 2 and on the right for T = 20 


algorithms for multiple subdomains is given by: for k = 1,2... and for 9 G (0,1] 
compute 

d t ul - d xx 
u\{x, 
u\ (X2 
u%{x 3 

and then for i = 2,1 

d t u% - d xx ii \ 
u\(x,$) 

u\(Xi-X,t) 

d x u$(xi,t) 

and finally for * = 4, 5 

d t u\ - d xx iii 

uH x ’ 0 ) 

d x Ui(xi-i,t) 

Ui(Xi,t) 

with g k (t) = g(0,t),g£(t) = g(5,t) for the first and last subdomains at the physical 
boundaries. The updated interface values for the next step are defined as 

9i(t) = @ u i(xi, t) + (1 - 6)g k i ~ 1 {t), i = 1 , 2 , 

9i(t) = 0u k i+1 {xi,t) + (1 - 9)g k ~ 1 (t), i = 3,4. 

We solve (2.1) using the above DNWR algorithm for different time windows for the 
same setting as in Al. Figure 2.5 gives the convergence curves for different values of 
the parameter 9 for T = 2 on the left, and T = 20 on the right. 

From the three numerical tests of the DNWR methods (Al, A2 and A3), it is ev¬ 
ident that the behavior of these algorithms are similar for smaller time windows. But 
we notice clearly faster convergence for the arrangement A3 for large time windows. 
We therefore focus on the third version (A3) of the DNWR algorithms, and formally 
define the DNWR method for the general parabolic model problem (1.1) for multiple 
subdomains in the next subsection. 


u k 

= 0, 

in 

D3, t > 0, 

0) 

= U 0 (x), 

in 



= gtW, 

t > 0, 


= g k 3 ~\t), 

t > 0, 


o, 


in flj, t > 0, 

= 

u 0 (x), 


in f 

= 

gtHt), 


t > 0, 

= 

9 x u i+l( x ii 

t), 

t > 0, 


0, 


in f2 i,t > 0, 

= 

uo{x), 


in f2j, 

= 

dx^i — 1 [Xi — 1 


t > 0, 

= 

gt\t), 


t > 0, 


( 2 . 2 ) 


(2.3) 


(2.4) 
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Fig. 2.5: Arrangement A3: Convergence of multi-subdomain DNWR for various re¬ 
laxation parameters 9, on the left for T = 2 and on the right for T = 20 


2.2. DNWR algorithm. We now formally define the Dirichlet-Neumann Wave¬ 
form Relaxation method for the model problem (1.1) on the space-time domain 
f l x (0, T) with Dirichlet data given on dfl. Suppose the spatial domain ff is par¬ 
titioned into 2m + 1 non-overlapping subdomains ffj, i = 1 ,..., 2m + 1 without any 
cross-points, as illustrated in Figure 2.6. We denote by u t the restriction of the so¬ 
lution u of (1.1) to fii. For i = 1,..., 2m, set T.; := df2j D dfl;+i- We further define 
To = F 2 m+i = 0- We denote by n l ,j the unit outward normal for f2,; on the interface 
rj, j = i — l,i (for f 2 i,fl 2 m+i we have only m : 2 and n 2 m+i, 2 m respectively). In 
Figure 2.6, D and N in red denote the Dirichlet and Neumann boundary conditions 
respectively along the interfaces as in the arrangement A3. 

The DNWR algorithm starts with initial Dirichlet traces g®(x,t) along the in¬ 
terfaces Tj x (0, T), i = 1,..., 2m, and then performs the following computation for 
k = 1,2,... 


d t u k m+1 - V • (/c(®,t)Vu^ +1 ) 

= fi 

in , 

U m+l( X i 0) 

= u 0 (x), 

in Q m _|_i, 

^m+1 

= 9i 

on dQ fl d£lrn+ 1 , 

4+1 

II 

<s> - ?s- 

1 

on I\, i = m, m + 1 


(2.5) 


and then for m > i > 1 and m + 2 < j < 2m + 1 


d t u^=V ■ (n(x,t)S7u’l) + /, in f 
u k {x, 0)=uo(x), inflj, 


u i = 9i \ 

d r ,, u k =-d, 




d t u k =V ■ (k(x, t)Vu^) +/, in Clj, 
u k (x,0)=u o {x) 1 in Qi 

on dfb \ r i; d„ ^u^-dn v!}^, 
k -1 


on r i 


rii,i +1 “i — ^rn+i.i “i+l) 

with the update conditions along the interfaces 


U j 9j > 


'3i 

onr.,_i, 
on df lj \ Fj_i, 
( 2 . 6 ) 


9i(x,t)=9u? |r iX (o,T) + (1 - 0)g i (x,t), 1 <i<m, 

9j( x i t) = Ou,j +1 |r iX (o,T) +(1 -0)g^~ 1 (x,t), m+l<j<2m, 


(2.7) 


where 9 G (0,1], and for i = 1,..., m and j = m + 2,..., 2m + 1 


~k -1 
9i = 


k -1 

. 9i -1 ! 


on df l (~l df li 
on Ti-i, 


9j 1 = 


on df2 (~l df2 j, 


g k j \ °n rj, 
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Fig. 2.6: Splitting into many non-overlapping strip-like subdomains 


Remark 2.1. The DNWR algorithm (2.5)-(2.6)-(2.7) is defined for an odd num¬ 
ber of subdomains. In case of an even number of subdomains 2m + 2, we treat in a 
similar way as above the first 2m +1 subdomains, keeping the last one aside. Then for 
the last subdomain, we apply a Neumann transmission condition along the interface 
r 2 m+i an d a Dirichlet boundary condition along the physical boundary. 


2.3. Convergence analysis. We present the convergence result of the DNWR 
algorithm (2.5)-(2.6)-(2.7) for the ID heat equation with n(x,t) = v. We split the 
domain Cl := (0, L) into non-overlapping subdomains Cli := {x%-\, Xi), i = 1,..., 2m+ 
1, and define the subdomain length hi := Xi—Xi-i. Also, define the physical boundary 
conditions as u(0,t) = go(t) and u(L,t) — gL(t), which in turn become zeros as 
we consider the error equations, f(x,t) = 0,go(t) = gL(t) = 0 = uq(x). We take 
{9°MZ " as initial guesses along the interfaces {x = Xi} x (0, T), and for sake of 
consistency we denote g^t) = g k m+ i = 0 for all k corresponding to the physical 
boundaries. Denoting by { w k (i)}) m | for k = 1,2... the Neumann traces along the 
interfaces, we compute 


“m+ 1 ^ 0 ) 
U m+li X mi t) 
u m+1 ( x m+ 1 1 1) 


0, X £ D m _)_i, 

0, X £ S2 m _|_i, 

9 k m ~ + \(t). 


( 2 . 8 ) 


and then for m > i > 1 and m + 2 < j < 2m + 1 


d t ui - vd xx u\ =0, ie!l„ 

u*(a;,0)=0, x€fli, 

u i ( x i-i, t) = g k I\ (t)> 
d x u'f{x i ,t) = w 1 f{t), 


d t u k - ud xx u k =0, x € Cl j, 

Uj(x, 0)=0, x £ Clj, 

u j(xj,t)=gj{t), 


(2.9) 


and finally the update conditions with the parameter 9 £ (0,1] 


g k {t) = 6u k i {x i ,t) + {l-0)g k 1 (t), w k {t) = d x u k +1 ( Xi ,t), l <i <m, 
g k (t) = 9u k +1 (xj,t) + (l-9)g k ~ 1 (t),w k (t) = -d x u k (x j ,t) 1 to + 1 < j < 2m. 

( 2 . 10 ) 
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We have the following main convergence result for DNWR for the heat equation and 
the proof will be given in Section 5. 

Theorem 2.2 (Convergence of DNWR for multiple subdomains). For 9 = 1/2 
and T > 0 fixed, the DNWR algorithm (2.8)-(2.9)-(2.10) for multiple subdomains 
of unequal sizes hi, 1 < i < 2m + 1 with /i max := maxi<j< 2 m+i hi and h min := 
mini<j< 2 m +i hi converges superlinearly with the estimate 


max 

l<i<2m 


9i IU”(o,T)< [2m -3 + 


2K 


t-ra +1 


erfc 


/ kh m i n \ 


max 


\ 2 \fvT) !<^< 2 m 


9? II 


L°°(0,T) 


Remark 2.3. The estimate in Theorem 2.2 holds for an odd number of subdo¬ 
mains. In case of an even number of subdomains 2m + 2, we define the algorithm as 
given in Remark 2.1, and the estimate will be of the form 


max || gf \\ l ™(o,t)< 2m - 1 + 


l<K2m+l 


2K 


t-m+1 


erfc ( max || || ioo(0 T s . 

\2\fvT) l<i<2m+l y 

( 2 . 11 ) 


2m — 1 appears in (2.11) as we calculate the estimate for 2m + 3 subdomains. 


3. DNWR for hyperbolic problems. In this section we define the Dirichlet- 
Neumann Waveform Relaxation method with many subdomains for the model prob¬ 
lem (1.2) on the space-time domain 12 x (0, T). This can be treated as a generalization 
of the DNWR algorithm for two subdomains, for which see [9]. As before, there are 
three possible arrangements: Al, A2 and A3. However, for consistency we formally 
define below the DNWR algorithm for the arrangement A3. 


3.1. DNWR algorithm. Suppose the spatial domain 12 is partitioned into 2m+ 
1 non-overlapping subdomains 1 h,i = 1 ,...,2 m + 1, without any cross-points as 
illustrated in Figure 2.6. For i = 1,..., 2m, set Tj := 912j fl 912j + i. For consistency, 
we set r 0 = r 2m+ i = 0. We denote by Ui the restriction of the solution u of (1.2) 
to 12 1 , and by n^j the unit outward normal for 12^ on the interface Tj,j = i — l,i 
(for 12i, 12 2m +i we have only m j2 and n 2m+ i i2m respectively). We define the DNWR 
method as in the arrangement A3, but one can also consider as in Al and A2. 

Given initial Dirichlet traces g°(x, t ) along the interfaces x (0, T), i = 1,..., 2m, 
the DNWR algorithm consists of the following computation for k = 1,2,... 


c 2 (x)Au^ n+1 

= /, 

in , 

i + i(*. 0 ) 

= v 0 [x), 

in , 

dtu k m+ 1 (*, 0 ) 

= w 0 {x), 

in , 

U m+ 1 

= 9 , 

on dfl fl d£lrn+ 1 , 

ut+i 

1 

^ -CO 

II 

on I\, i = m, m + 1 


and then for m > i > 1 and m + 2 < j < 2m + 1 


d tt u , f = c 2 {x)Au 1 l + /, 

in 12 j, 

d tt u k j =c i {x)Au k j + /, 

in 12 j, 

«i(*, 0 )=i>o(®), 

in 12 j, 

Uj{x, 0) = vo(x), 
d t uf(x, 0 ) = w o (x), 

in 12 j, 

d t u*(x,0)=w 0 (x), 

in 12 j, 

in 12 j, 

1 

« •«» 

II 

-Sfi-ea 

3 

011 <912,; \ Fj, 

d njJ -i u j = —dn J _ 1J Uj_ i , 

on r^i, 

drii' i+1 uf — d ni+l i U i+l , 

on r, : , 

uk j=9 k 3 ~\ 

on 912 j \ r 


(3.2) 
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with the update conditions along the interfaces 

gf(x,t) = 9u* |r iX (o,T) + (1 - 9)g k ~ l (x,t), 1 < i < m, 

9 k (x,t) = 9u k + i|r\,x(o,T) + (! - 0)g k ~ 1 (x,t), m+l<j<2m, 

where 9 G (0,1], and for i = 1,..., m and j = m + 2,..., 2m + 1 


“fc-i 


9i 


id-1 ’ 


on dfl D <911;, 
on IVi, 


-fc- 1 



on dfl (~l dflj , 
on Tj , 


3.2. Convergence analysis. We analyze the DNWR algorithm (3.1)-(3.2)-(3.3) 
for the ID wave equation with a constant wave speed, c(x) = c. Consider a splitting 
of the domain Q := (0, L) into non-overlapping subdomains 12; := i = 

1,..., 2m + 1, and define the subdomain length gt := Xi — Xi—\. As we consider the 
error equations, the physical boundary conditions u(0,t) = g$(t) and u(L,t) = PlA) 
become zeros along with /( x, t) = 0, vo(x) = 0 = wq(x). We take {<7°(i)} i=1 as initial 
guesses along the interfaces {x = Xi} x (0 ,T), and for sake of consistency we denote 
9o(t) = 92m+l = 0 for all k corresponding to the physical boundaries. Denoting by 
{w k (t)}^ for k = 1,2 ... the Neumann traces along the interfaces, we compute 


d ttUm+i - c 2 ®* 


d t u. 


l m+l 


ra +1 


■U k m+1 

= o, 

X G 


(x, 0) 

= 0 , 

X G 


(x, 0) 

= 0 , 

X G 

(3.4) 


/t 'm +1 


(^m 11') — 


k~i 

i/ra 

9m+l 


U m+1 {Xm+1 j t) 

and then for m > i > 1 and m + 2 < j < 2m + 1 


0 ), 

(*)■ 


d tt u k - c 2 d xx u k = 0, 
u?(ar,0) = 0, 
d t u k (x, 0) =0, 


G 12;, dttu k - c 2 d xx u k = 0, 


X G 12;, 
X G 12;, 


u'l{x i - l ,t)=g i _ 1 (t), 

,k (™ j.\ n .,k ( 


-d a 


u k Ax, 0) = 0, 
d t vh(x, 0) =0, 


X € 1 Ij, 

X € 12 j, 
X € 12 j, 


'j 


U k (x 


J-!,", ~j-J 

3 ,t)=g k ~ 1 {t), 


d x u?(xi,t) = wf(t), 

and finally the update conditions with the parameter 0 G (0,1] 


k -1 


g k (t) = 9u k (xi,t) + (l-9)g, 

g k (t) = 9u k +1 (xj,t) + (l-9)g k ~ 1 (t),Wj(t) = -d x u k (x j ,t), m+l<j< 2m. 


it) 

k -1 


(t) = d x u k +1 (x.i,t), 1 <i< 

k 


(3-5) 


(3.6) 


We now state the main convergence result for DNWR for the wave equation. The 
proof of Theorem 3.1 will also be given in Section 5. 

Theorem 3.1 (Convergence of DNWR for multiple subdomains). Let 9 = 1/2. 
Then the DNWR algorithm (3.4)-(3.5)-(3.6) converges in at most k + 1 iterations for 
multiple subdomains, if the time window length T satisfies T/k < h. m c, where c is 
the wave speed. 

4. Auxiliary results. We need a few auxiliary results related to Laplace trans¬ 
form to prove our main convergence results. We define the convolution of two functions 
g , w : (0, oo) —> R by 


{g*w)(t):= g(t - T)w(r)dT. 

J o 
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Lemma 4.1. Let g and w be two real-valued functions in (0, oo) with w(s) = 
C{w{t)} the Laplace transform of w. Then for t £ (0 ,T), we have the following 
properties: 

1. If g(t) > 0 and w(t) > 0, then (g * w)[t) > 0. 

US' * HUro.t) ^ II5||l 1 (o,t)IMIl 1 (o,t)- 

3 . |(5*w)(<)|< ||5lU- ( o,T)/ 0 T |w(r)|dr. 

4■ fg w(r)dr = ( H*w){t ) = Cr x H(t) being the Heaviside step function. 

Proof. The proofs follow directly from the definitions. □ 

Lemma 4.2 (Limit). Let, w(t) be a continuous and L 1 -integrable function on 
(0, oo) with w(t) > 0 for all t > 0, and w(s) = £{w;(f)} be its Laplace transform. 
Then, for r > 0, we have the bound 



lim w(s). 

s—>• 0 + 


Proof. For a proof of this lemma, see [10]. □ 

Lemma 4.3 (Positivity). Let /3 > a > 0 and s be a complex variable. Then, for 
t e ( 0 , oo) 


v(t) ■= £ 1 


sinh(gyi) ) > 

sinh(/3i/s) J _ 


and ip(t) := C 1 


cosh(aVi) \ > 
cosh(/3-/s) J _ 


Proof. For a proof, we again cite [10]. □ 

Lemma 4.4. Let a,/3> 0 be two real numbers and s be a complex variable. Set 


_ = siiih((q - P)yfl) 

^ cosli(a- v /s) sinh(/3- v /s) 


(4.1) 


Then 



£ 1 {*(•)} (dl* - 


lim x(s) = 

s—>- 0 + 


a — /3 


Proof. There are two possibilities: a > /3 or a < /3. In either case, we need 
Lemma 4.3 to prove positivity of the expression. For a detailed proof, see [10]. □ 

5. Proof of main results. We now prove the main convergence results for the 
DNWR algorithm applied to the heat and wave equations, stated in Section 2 and 3. 

5.1. Proof of Theorem 2.2. We start by applying the Laplace transform to 
the homogeneous Dirichlet subproblem in (2.8), and obtain 

s Wrc+l — = 0) Wrc+l ( x m, s) = CJ m (s), W m +1 ( x m+l , s) = 5m+!( S ) J 


Defining cq := sinh [hi^/s/is'j and ji := cosh | 'hiyJs/i^j , the subproblem (2.8) solu¬ 
tion becomes 

u^ +1 (x, s)=-^^ (j)5n+i ( s ) sinh {ix-x m )y/sjvj sinh ({x m+1 -x)y/sjvf) ' 
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Similarly the solutions of the subproblems (2.9) in Laplace space are 

Ui(x, s) = ~ Xi-^s/sjv) + ^-9 k i-1 cosh((xj - x)sfsjv), 

Uj{x, s ) = ^ _1 cosh((a; - Xj-Jy/s/v) + sinh((a; i - x)^/sjv). 


for 1 < * < m and m + 2 < j < 2m + 1. Therefore for 9 = 1/2 the update conditions 
( 2 . 10 ) become 


9 k i = 


1 n ^ 1 + In fc - 1 + 

9i -i ' 2 9i ^ 


27, 


^ v/^’ 


9 k 


g .7 + l I 1 7 fc-l 

2 7i+l yftjv 


i= (5- 1 ) 


for l< 2 <m, m + 1 < j < 2m and 


7i + l — — ' " L V V CTm+ 1^" L ' O 

&n+: i; rn + 2 < j < 2m. 

(5-2) 

We choose g* := 7 ^, wf := -^=< 7 *, 1 < i < rn and g* := 'y j+1 g^,w^ := -^L=a j+ll m+ 
1 < j < 2 in with 70 = 72 m +2 = 1 in (5.1)-(5.2) to get 


,fc__ ys~ < T i+i afc-i 
* V ^ 7*+i 

,t,fe _ VfZE^-l _ ^~ 7m + l flfc -1 . 

'"m+l a m+1 9m V 1 / CT m+1 9ra+l 1 


1 «j? +1 ,l < * < m- 1 ; 


rf.k _ _ /JTm+i "fc-1 
m V cr m+ i Vm 


y s/v ~k- 1 
fm + 1 “ m + t’ 


/«7^ 


+ 3 % + 3 ™?, 1 < * < m; ^= 3^1 + 14 + 2 ^ 4+1 ’ m + 1 < j < 2 m, 


and for 1 < z < m — 1 and m + 2 < j < 2 m, 




7 T,fc . .,7,fc — <*mltm+l Tzk— 1 


7«7i+l 1 tTi+17i+l *+l ’ m 7m°'m+l i ' m ' 0’m + l7m+2 

„T.k _ °~m+2 7:fc —1 _ 7m+lfr’ m +2 . n -7,k _ a j+l „T.k _ a j a j+ 1 7 ^k—l 

777-*-! - - , 9m cr m ^-i7m+2 ! 7 'VAfTA 7 — 1 i i 27 


- fc -1 

Sm+ 1 ' 


m + 1 1m7m+l' 


7 3 -o-i 


7j7j + l 


Therefore in matrix form, the system becomes 


( ~9 f \ 


A 


„T,k 


®i+i 

ffm+1 


= B 




9 k m 1 

Sm+i 


, where A 


1 0 

0 Lji 


4m X 4m 
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K, = 


<? i 


272---7mO'm + l7m+2 
<7 2 


K a = 


273 -7m1m+l7m+2 

(7m — 1 

2 7 m 

CTm + l7m+2 


(7m 

2cr 

m + l7m+2 

£7^ 

i + l,m + 2 

2(7 „ 

i+l7m+2 

(7 j 

7i+37m + l 


— C r 2m7m + 1 

- 


,tf 3 = 


C m + 2 
27m<7m+l 

_°~m-|-3_ 

27mO'm + l7m+2 


27m "m+1 7m+2 ■ ■ -72m-1 

o - CT2m+1 - 

27 m(Tm+ 17 m+ 2 ..- 72 m . 


2 7m+3 
7m + 2,m + 3 
27m + 27m+3 


-^ 2 mffm +2 
V3- 


27 m +4 


2o- m + l7m+2'>'m+3---72m-l 27m-|-27m+ 3 7m+4 ■ ■ -72m-1 
_ — "2m+l7m+l _ _ — "2m+l "m + 2 _ 


72m — 1,2m 1 

272m —172 m 27 2 m-f-l 

_ "2m-l"2m+l 72m,2m+l 

272m -172m 27 2m 72m+l 


2lTm+l7m + 27m+3---72m 27 m+ 2 7m+ 3 7m+4 ■ ■-72m 

with 7 ^ = cosh ((h^ — hj)\Js/v'j , <Tjj = sinh ((h^ — hj)\Js/v'j . Therefore we can 
write the first equation as: 


ai 7172 717273 a 1 


_ CTlCTm *fc~ 1 J l7m-fl ~fc—1 

7l72---7m-l7m "m-l 7l...7m(Tm+l 


ff i a fc -i 

7 i...7 m o- m +i ^rn+1* 


(5.4) 


So with the substitution v * := 2 fc q k g^, 7 = cosh ( h m i n y/s/is ), we get 


£.fc _ 771,2 gfc-l 70-10-3 7 .fc-l_ 7 <llJm +.fc~l 7 g~l 7 m+l £.fc-l | 70-1 £- fc ~ 1 

1 7172 1 717273 2 7l72-.-7rn-l7m "1-1 7l...7mffm + l m T 71 ■ ■ -7m"m + 1 "1+ 1 


m +1 




• k-l 


i=l 


Now using Lemma 4.3 and 4.2, we obtain 

700 cosh(/i m i n \f~sjv) cosh I (hi — h,2)\Js/v\ 

/ \dn(t)\dt < lim - -j= - -== - < 1 , 

Jo ’ s-s-o+ cosli(/i2 \/s/v) cosh(hi^/s/v) 


so that using Lemma 4.1, part 3, we get 

Jr -1 (d^vt 1 ^)) 

Moreover, we write for j = 2,..., m — 1 

7 


L°° ( 0 ,T) 


< IK*" 1 ! 


L°° ( 0 ,T) ■ 


rf l,J = 


x _ Q'lg'j+i _ 1 } = 

72---T j \ 7i7l+i 


\ 7 cosh ((hi — hj+ijy/s/v 

) 72 • ■ ■ 717171+1 


7 


72 • • ■ 7l 


cosh ((h m in + hi — hj+i)y/s/v^ + cosh ((h min - hi + hj + i)y/sj 

272 ■ 7l7l7l+i 


7 


72■ ■ • 71 


Now —hj+i < hmin + hi — hj+i = hi + h m i n — hj+i < hi, so for the first term, we 
choose the pairing 


1 


1 


cosh ((hmin + hi - hj+l)s/sjv^j 


72 • • • 7 j 2 cosh(hy + i \fsjv) 


cosh(hi-\/s/V) 
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and 


cosh ({h mia -hx+ h j+1 )y/aju\ 


72 • • • 7 j 2 cosh(/ii y/s/v) 
We therefore get by Lemma 4.3 and 4.2 


cosh(/ij+i y/s/v) 

r°° i , . m „ „ ,. (cosh((h mitl +h 1 -h j + 1 )y/s/v)+cosh[(h Ini n-h 1 +hj + 1 )y/s/v) \ 

In \d\ At)\at< lim ---;;------ H- 1 — 

Jo I UJWI — s ^o+ \ 2 72 —7i7l7j+i 72■■■73 y 


< 2 , 


so that 


C~ 


(^«f' 1 w)|| Loo( 0 ir) < 2 ||t;: 


fe —11 




Finally we write di m = -J— (1 - ffl7m+1 - l) = -2— . A+s)-2—, where 

J 72---7m V 7l°’m+l ) 72---7m V ' 72-”7m 


Ad(s) = - 


smh((/ii — h m +i)-^s7 


sinh(/i m _j_i s/y) cosh(^i -y/s/i') 5 


and write 


di 


1 


,m+l — 


/ ovy _ 
ff m +l 72 ■ • - 7m V71CT 


1 + 1 = 


1 -A+s) 


^m+l T2 ■ • • Tm ^Vn+l T2 • • • Km 

( !——\ sinhf (hi — /imin)y s/v J 

hminy/s/v), K 2 (s) = - v , - 4^^. Note that both 

/ sinh(/iminV s/i/) cosh(/iiY's/i/) 

Ad(s) and K 2 (s) are of the form (4.1). So by Lemma 4.4, we obtain the bounds 

I hi | 


nOO 

/ \dl,m(t)\ 

JO 


dt. < lim 


7 


s ^o+ 72 ... 7 „ 


(A' 1 (s) +1) < 


*'771+1 


and 
(*00 


roo 

/ |di,m+i(i)|<ft < lim 

Jo S-30H 


(7 1 


-S-0+ CT m+ l 72 . . . 7n 


(A+s) + l)< 


^min (hi h 


ra+1 \ ^min 


+ 1 = 




*'771+1 


We therefore get 


KOI 


L°°(0,T) 


< l + 2(m-2) + 


2 /in 


*-777+1 


■) i<Kfi 11 


< 


2 m-3+^^^ max || I|l«-(o,t) • 

ftm+l / l<J<2m 


We also get similar relations for other equations. By induction, we therefore get for 
all i 


KOI 


L°°( 0 ,T) 


< 


< ■■■ < 


2m-3+^^] max \\ v* +) |U°°(o,r) 
n-m+l J l<J< 2 m 

9 h \ k 

2m - 3 + max || 5 ?(-) ||l-(o,t) ■ (5-5) 

rtm+l / l<J<2m 


Setting f k (t) :=C 1 |^} we have 

9i(t) = ^ I (/ fc * v i) k = 2 i¥ / / fc (* ~ T ) v i( T ) dr , 
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from which it follows, using part 3 of Lemma 4.1 that 

llfl?(-)IU«(o,n < £ \f k (r)\dT. (5.6) 

By Lemma 4.3, f k (t) > 0 for all t. To obtain a bound for f k (r) dr , we first show 
that the function r k (t) = C~ 1 { 2 fe e _fc,lmin V^7^) is greater than or equal to f k (t) for all 
t > 0, and then bound j f J r;.(r) dr instead. Indeed, we have 


C {r k (t ) - f k (t)} = 2 k e - kh - 




s/v _|_ g ^min^/ s/l^^k 


= E (^e-^^sech 

l=i 


Also, by [1 



j h min ~-J 2 hl,„/vt 

Vnut 3 


(5.7) 


is a positive function for j = 1,..., fc. Thus, £ 1 sech fe (/t m in-\/E^')) i s 

a convolution of positive functions, and hence positive by part 1 of Lemma 4.1. This 
implies r k (t) — f k {t ) > 0 , so we deduce that 



dr < 



dr = £ 1 



(5 ' 8) 


where we expressed the second integral as an inverse Laplace transform using Lemma 
4.1, part 4, which we then evaluated using the following identity from [19]: 

£_1 (i'~ AV5 ) - ertc (iyi) • A >°- < 5 - 9 > 

Finally, we combine the above bound (5.8) with (5.5) and (5.6) to conclude the proof 
of the theorem. 

Remark 5.1. For equal subdomains, with subdomain size h, we can deduce a 
better estimate: 


max 

l<t< 2 m 


9 i 


IU°°(o,r)< (min {2m- 1, Q(h, v, T)}) k erfc 



max 

l<i<2 m 


II 9i I|l~(o,t), 


where Q{h, v, T) := 2 erfc ) + E“o 2 i+1 erfc • 

Here is the outline of the proof: For equal-length subdomains we have hi = h for 
i = 1 ,. .., 2 ro + 1 , so that = 7 , 0 , = a and 7 ^ = 1 , <Jij = 0 for all i,j in the matrix 
P in (5.3). Therefore, the first update equation (5.4) becomes 


9$ = 


2 7 


a k ~F 

2 in 


— n k ~F 

3 52 


27 


27 


m Um-l 


-fc -1 


27' 


771—1 *■ 


+ m — 2; (5.10) 
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which in turn becomes after setting v k (s) = 2 k j k g k (s), 


1 


v k = -v k ~ l - ( 1 - ^ 


7 

ra +1 

i =1 


I "2 


v!;- 1 - 


l 


v m—3 


l - 


l 


r 


~k -1 


1 


,2 j U m -1 2 U ™ 


c.k -1 


1 ^ 771+1 


(5.11) 


iVoie that the estimate from Theorem 2.2 also holds in this case with h max = h m +±: 

max || g\ |U~(o ) t ) < (2m-l) fc erfc (-^=) max || 5 ° I|l°°(o,t) ■ (5-12) 

We now present a different estimate starting from (5.11). By back transfoi'ming 
into the time domain we obtain v k (t) = ( r * * 0 )- Now part 3 of Lemma 

f.l yields 

IK fc IU-(0,n ^ E ll«? _1 IU“( 0 ,T) [ T \n{r)\dT 

i=1 J o 

II V k ~ 1 ||l°°(0,T) E [ \ r i( T )\ dT - ( 5 - 13 ) 

i=l J ° 


< max 
l<£<m+l 


Set f k (t) := £ 1 . Since lim s _>o+ 11 — 1 <2, we have from (5.13) 

TTTL~\~ 1 r Jf 1 r Jf 

[ \ r i{ T )\dr < f (fl(t) + 2 + 2fl(t) + ■ ■ ■ + 2f m _s(f) + fm- 2 {t) + fm-lit)) dt. 
i=1 J ° J ° 


Therefore using the inequality (5.8), / Q T fk(t)dt < 2 fc erfc |' —^=? J, the above expression 
is bounded by 


m+1 r j 

T-.L 


|ri(r)|dr < 2 erfc 


2\fuT J 


) + E 2l+1 erfc ^ Q( h > v > T )> 


i =0 


where Q{h,v : T) := 2 erfc ^= j + 2 t+1 erfc ^ .,^‘ /r ) ■ So we get the second 
estimate as 


max II sf IIl»(o,t)< erfc 

l<i<2m 



max || g° ||l°=(o,t) . 


(5.14) 


The result follows combining the two estimates (5.12) and (5.14). 

We compare the two estimates (5.12) and (5.14) in Figure 5.1 for u = 1. The 
region below the red curve is where the estimate (5.12) is more accurate than (5.14). 

5.2. Proof of Theorem 3.1. We start by applying the Laplace transform to 
the homogeneous Dirichlet subproblem in (3.4) to get 




L m+1 


~c 2 u k 


m-\-l,xx 


= 0 , 


i 'ra +1 


(x m ,s) — g m 00 , u m+1 (x m +i, s) — g m +i{s) 


Define p, := sinh (his/c) and A i := cosh (his/c). Then the subproblem (3.4) solution 
becomes 

«m+iO> s ) = — — {9^n+i(s) sinh ((a: - x m )s/c ) + ^ _ 1 (s) sinh {(x m+1 - x)s/c)) . 
Pm +1 
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| Q(h,T)=2m-1 I I Q < 2m-1 | 



Fig. 5.1: Comparison of the two estimates (5.12) and (5.14) for v = 1 


The solutions of the subproblems (3.5) in Laplace space are 


A i s/c 

u k Ax , s)=j-g k ~ 1 cosh((a; - Xj-i)s/c ) + 


u k (x, sinh((:r — a: i _ 1 )s/c) + cosh((xi — x)s/c ), 1 < i < m, 

’ i 

A 3 - s/c 


sinh((xj — x)s/c), m + 2 < j < 2m + 1 . 


Therefore for 0 = 1/2 the update conditions (3.6) become 


3 ? = hat 1 + 1 < * < 

9j = + tit 1 + 217773-+i, m + 1 < 


; < 2m. 


(5.15) 


and 


= _*£<±i ^- 1 + i , 1 < i < TO - 1 ; wt = 

1 C Ai +1 — — ’ m C Pm.- H 0m- 1-1 5 

— S / C 1 _ s A m+ 1 -fc-1 . 7 7,/c _ J_ 7/r _ 

^m+i p m+ i 9m c pm+i 9m-\-l i W j A 7 W j-1 c A, 


CPm+r m Pm + 1 * 

'’Ll - § Lg? -1 , m + 2 < j < 2m, 


(5.16) 

We choose := \g k ,w k := 1 < i < m and g* := A ; . ;</*.//;? := ^ pj+i,m + 

1 < j < 2m with Ao = A 2 m +2 = 1 hr the corresponding equations of (5.15)-(5.16) to 
get 


50 ? 1 + 5 W ?> 1 <i<m; 


0 ? = 2ATT3--1 1 + 1 + 

3? = H* + \g k ~ l + 5x7773?+i > ™ + l<j< 2m, 


and for 1 < i < m — 1 and m + 2 < j < 2 m, 


,7,fc PiPi + l -k— 1 J_ Pi -k . jk PmKn + i_-k — l j_Pm_A”! 

Wi AiA i+ i 3i + W m XmPm + i9m p m+1 A m+2 3m+l> 

-fe _ Pm + 2 -fc-1 __ A m + ip m + 2 -fc-1 . -k _ Pj + 1 jfc _ PjPj + l ^fc-l 

m +! A m p m+ ii/rn Pm+ 1 A m +2 i/m+1) w j \ jPj w j-l XjXj+i 9j 
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Therefore we can write the system in matrix form as in the Heat equation case 


(\ 



—k 


— k—1 

9m 

—k 

= P 

9m 

-k-1 

9m -\-1 


9m +1 

V 92m ) 


\92m) 


where P is in the same form as in (5.3) with 7 ij and ovy replaced by A ij := 
cosh ((hi — hj)s/c) , pij := sinh ((hi — hj)s/c) respectively. Therefore the updating 
conditions become 


ra+1 


0 +1 


9i( s ) = k, l 9i X (s), 1 <i<m\ 9ji s ) = 1 (s), m + 1 < j < 2m, 


l=i— 1 


l=m 


where fci >0 = /c 2 m, 2 m+i = 0 , and for i + 1 < l < m k^i-i = 7 ^-, = 


pipi+i 


ki.m, — 


Pi X m+l t. _ Pi 

2Xi...\ mPm+1 2A,...A m p m + 1 


(5.17) 

A c»+i k. , _ 
2A,A i+1 > n 'V 


2AiA <+ i...A,+i’ n '*- m 

TO + 1 < / jg j - 1 kjj = 2AJa7+i ’ = ~ 2Xi\?+i. P 'x j+ i ’ = 2 A J 1 +1 ’ %> m +! = 

Pj+i Am + i t _ _ Pj + i 

1 •v'l.m. — 


for 1 < i < to, and for 

l 


-75 -r-r- ,n>im — 75 - £l+l -- f or m + 1 < j < 2TO. AlsO, k m m - 1 = 

2p m + lA m + 2...Aj + l ’ j’ ln 2p m +lAm + 2...A, + l J — ’ "‘.111 

1 7 — gi2±uia £ ,. — Pm. „ nr i £. 

2A m ’ A ' m 5 m — 2p m _(_iA m ’ n '7n,m+l — 2p m _i_iA m A, m+l,m — 


Pm -\-2 


2p m + lA m + 2 ’ ^ m + 1 . m + 1 _ 

2 p m +iT + + 2 ’ ^m+i,m +2 = 2A 1 +2 ■ So by induction on (5.17) we can write for 1 < i < 2m 


2m 

9i(s) = ^2ptj (kl,l,h,2, ■ ■ ■ ,k2m,2rn-lik 2 m,2m) 5* _ "(s), (5.18) 

j =1 

where the coefficients p r - j are either zero or homogeneous polynomials of degree n. 
Now expanding hyperbolic functions into infinite binomial series, we obtain for i +1 < 
l < m and 1 < i < m 


ki i 


cosh ((hj - h l+ i)s/c) 

2 cosh (his/c) cosh(/ii+i s/c] 


— I e -2his/c _|_ e ~2h i+1 s/c 


1=1 


(-l)'e- 


-y^(_l)" e -2fti+in S /c 


oo oo 

^ ^ 2(/fei+n/ij+i)s/c 

/=ln=l 


ki,i = - 


sinh ( his/c ) sinh (hi+is/c) 


2 cosh (his/c) cosh(/i i+ i s/c) ... cosh(/q + is/c) 

e ~2his/c _ e ~2hi +1 s/c e -2(/j i +/j i+ i)s/c^ ^ _|_ g -2/i„s/c^ 




k% i—l — 


2 cosh (his/c) 


— e ~his/c 


1 + E<-!) 


I ^—2hils/c 


1=1 
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sinh (his/c) cosh (h m +is/c) 


2 cosh (his/c).. .cosh (h m s/c) sinh (hm+is/c) 

_ _2 m— i ~ 1 e~ ( hi + 1 ^ - \-h m )s/c ^ _ g-2 his/c _|_ e ~2h m+1 s/c 


_g — 2(hi+h m + 


- 2 h m+1 8 /c\ Y[(l + e- 2h,s/c ^ , 

l=i 


sinh (hi s/c) 

i,m+i 2 cosh(/ijs/c)... cosh (h m s/c) sinh (h m+ is/c) 

= 2 m - i e~( hi+1 ~ { - \-h m+ i)s/c ^ _ g-2 his/c^j ^ _ e ~2h m+1 s/c^ 1 ^ _|_ e ~ 2h ‘ s / c ^ 


-1 


l=i 


kra.m — 


sinh {(h m+ i - h m )s/c) 

2 cosh (hms/c) sinh(/i m +is/c) 


(' 


_ , g-2 h m s/c _ e -2h m+ 


is/c ) [1 


+ J2(-l) l e~' 2hmls/c + J2e~ 2hm+ins/c + 


/g-2(//i m +n/i m .|-i)s/c 


1=1 


n= 1 


1=1n=l 


k 


m,m-\-l 


sinh (hms/c) 

2 cosh (h m s/c) sinh(/i m+ i s/c) 


^g hjn-\-\s/c 


g—(2h m +h m _j_i 



[i 


oo oo oo oo 

+£(-1 ) l e- 2h - u ' c _|_ ^g-2/i m +ins/c _|_ EE j_|jlg-2(l/i m +nfi m+ i)s/c 

/—1 n=l 1=1n=l 

The argument also holds similarly for other terms. Now using these expressions we 
can write (5.18) as 


gHs) = 


2m 

E ffc 

1=1 


* j 


Wa?( 


(5.19) 


where -(s) are linear combinations of terms of the form e sz with z > khi/c for 
some l G {1,2,... 2m +1}. We now recall the shifting property of Laplace transform: 

C- 1 {e- QS /(s)} = H[t - a) fit - a), (5.20) 


where Hit) := < ’ ^ > is the Heaviside step function. We use (5.20) to back 

\o, t<0. 


transform (5.19) and obtain 

9 ?(t)=g°(t- 


khi 

c 


H 


khi i 

t -— ) + other terms, 


for some j G {1,2,... 2m} and l G {1,2,... 2m + 1}. Thus for T < fc/i m ; n /c, we get 
9i(t) = 0 for all i, and the conclusion follows. 

Remark 5.2. The shifting property of Laplace transfoi'm (5.20) is the reason 
behind the finite step convergence of the DNWR for 9 = 1/2. The right hand side of 
(5.20) becomes identically zero for t < a, so that for sufficiently small time window 
length T (e.g., T < a) the error becomes zero and leads to convergence in the next 
iteration. In Figure 5.2 we plot Tr x |/(s) j w /(£) = sin (t) on the left, and show 
the effect of time-shifting on the right. 
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Fig. 5.2: Example of time-shifting for the function /(f) = sin(f): C 1 |/(s) j 011 the 
left, and £ _1 je _s /(s)j on the right 


6. Analysis of DNWR algorithm for Wave equation in 2D. We now 

formulate and analyze the DNWR algorithm for the two-dimensional wave equation 

d tt u — c 2 Au = /(x,y,t), ( x,y ) G ft = (l,L) x (0, tt), t G (0,T] 


with initial condition u( x, y , 0) = vo(x, y), dtu(x, y, 0) = Wo(x , y) and Dirichlet bound¬ 
ary conditions u(x,y,t) = g(x,y,t),(x,y) G <9f2. To dehne the DNWR algorithm, 
we decompose ft into strips of the form f2j = (Xj_i,Xj) x (0,7r), l = Xq < X\ < 
••• < X 2 m+i = L. We define the subdomain width hi := X; — Xi-i, and fi m in := 
mini<i< 2 m +i hi. Also we directly consider the error equations with /(x, y,t) = 
0,i>o(x,y) = 0 = wo(x,y) and homogeneous Dirichlet boundary conditions. Given 
initial guesses {y°(y,t)}(_ 1 along the interface {x = Xj}, the DNWR algorithm, as a 
particular case of (3.1)-(3.2)-(3.3), is given by performing iteratively for k = 1,2,... 


4+1 - c 2 A?4 +1 

= o, 

in 


= o, 

in 

9 t M^+i(x,2/,0) 

= o, 

in 

M m+1 { x m i 2/ 7 

= 



= g^+iiyA), 


“m+ifoM) 

o' 

II 

£ 

+ 

£ 

s 

II 



and then denoting by {w k (y, t)}^™ the Neumann traces along the interfaces, calculate 
for m > i > 1 and m + 2 < j < 2 m + 1 


9ttU?-c?Auf = 0, 

in ft. 

(^7 2/ 5 0) = 0, 

in ftj 

d t u k (x,y, 0)=0, 

in ftj 

u k i {x i - ll y,t) = g k ll{y,t), 


d x u k (x i ,y,t)=w k (y 1 t), 



u( C (x,0,t)=wf(x, 7T,t) = 0, 


dttUj - c 2 A u k = 0, in ft.,-, 

u k Ax,y, 0) = 0, in ft,-, 

d t Uj(x, y, 0) = 0, in ft,-, 

d x u k (x,_i, y, t) = w k _ x {y, t), 
u k (xj,y,t)=gj{y,t), 
u k (x , 0, t) =u k (x, 7r, t) = 0, 

( 6 . 2 ) 


with the update conditions for 1 < i < m and m + 1 < j < 2m 


9i(y,t) = Qu\ (xi,y,t) + (1 — 0)g k 1 (y, t), w k (y, t) = d x u k +1 (, x u y,t), 
9j(y, i) = 0u k +1 (xj,y, t) + (1 - 0)g k ~ 1 (y, t),w k (y, t) = -d x u k (xj,y, t), 


(6.3) 
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where 9 € (0,1]. 

We perform a Fourier transform along the y direction to reduce the original prob¬ 
lem into a collection of one-dimensional problems. Using a Fourier sine series along 
the y-direction, we get 


y> t ) = '52 u i( x ’ n > ^ sin ( n y ) 


n> 1 


where 


2 f n 

Ui (x , n,t) = — v!f (x, r), t) sin (nrj)di). 

n Jo 


The equation (6.1) therefore becomes a sequence of ID equations for each n, 


d 2 Uj+ 1, 

dt 2 


(x, n, t) - c — ^ (x, n, f) + c n U„ +1 (a;, n, t) = 0 

with the boundary conditions for n, t). We now define 

X(a, /?, f) := U” 1 1exp (-/3Vs 2 + a 2 ) } , Re(s) > 0, 


(6.4) 


(6.5) 


with s being the Laplace variable. Before presenting the main convergence theorem, 
we prove the following auxiliary result. 

Lemma 6.1. We have the identity: 


X(a,/M) = 


t>0, 

0 < t < p, 


0 


where S is the dirac delta function and J\ is the Bessel function of first order given 
by 

i r 

J\(z) = — / cos (z sin <p — ip) dip. 


Proof. Using the change of variable r = \Js 2 + a 2 we write 
e~ pr = e _/3s - (e _/3s - e~ Pr ). 

From the table [21, p. 245] we get 

Also on page 263 of [21] we find 


C 


-l |g—/3s _ e -/ 3 r| = ( 02^2 Jl (“V^ ~ff i ) > < > Z 3 ’ 

0 < t < (3. 


( 6 . 6 ) 


(6.7) 


Subtracting (6.7) from (6.6) we obtain the expected inverse Laplace transform. □ 
Now we are ready to prove the convergence result for DNWR in 2D: 
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Theorem 6.2 (Convergence of DNWR in 2D). Let 6 = 1/2. For T > 0 fixed, the 
DNWR algorithm (6.1)-(6.2)-(6.3) converges in k + 1 iterations, if the time window 
length T satisfies T/k < h m - ln /c, c being the wave speed. 

Proof. We take Laplace transforms in t of (6.4) to get 


(^cV)C-c' 


d 2 U: 


ra +1 


dx 2 


= 0 . 


We obtain similar sequence of equations by applying Fourier series first and then 
Laplace transform to (6.2). We now treat each n as in the one-dimensional analysis 
in the proof of Theorem 3.1, where the recurrence relation (5.19) of the form 


9i ( s ) = ( s ) 9j («) (6-8) 

3 

become for each n = 1,2,... 

6?(n, «)=£ (y S 2 + c 2 n 2 ) G°(n, a). (6.9) 

3 


In the equation (6.8) f// (s) are linear combination of terms of the form e~ es for 
q > khi/c for some l £ { 1,2,... 2 m + 1}. Therefore the coefficients f * ■ (\/s 2 + c 2 n 2 ) 
are sum of exponential functions of the form e ~eVs 2 +c 2 n 2 f or ^ Hence we use 

the definition of x i n (6-5) to take the inverse Laplace transform of (6.9), and obtain 


Gi(n,t ) = 

3 l 

with qij > kh m [ n /c. So it is straightforward that for t < kh m - m /c, G 1 * (n, t) = 0 for each 
n, since the function \ is zero there by Lemma 6.1. Therefore the update functions 
g^{y, t), given by g*(y, t) = XZ n >i G\(n, t) sin (ny) are also zero for all i £ { 1 ,..., 2m}. 
Hence one more iteration produces the desired solution on the entire domain. □ 

7. Numerical Experiments. We show some experiments for the DNWR al¬ 
gorithm in the spatial domain 12 = (0,5), for the problem dtu = Am, with initial 
condition uo(x) = x(5 — x ) and boundary conditions u(Q,t) = t 2 ,u(5,t) = te - *. We 
discretize the heat equation using standard centered finite differences in space and 
backward Euler in time on a grid with A* = 2 x 10~ 2 and At = 4 x 10 -3 . In the 
first experiment we apply the DNWR for a decomposition into five subdomains and 
for three different time windows T = 0.2, T = 2 and T = 8, whereas for a fixed time 
T = 2 we run another experiment for three to six equal subdomains. In Figure 7.1, 
on the left panel, we show the convergence estimates in the five-subdomain case as a 
function of time T, whereas on the right panel, we show the convergence for T = 2 as 
we vary the number of subdomains. We observe superlinear convergence as predicted 
by Theorem 2.2, and for small T the estimate is quite sharp. We also see that the 
convergence slows down as the number of subdomains is increased. 

We now consider the following model wave equation to see the convergence be¬ 
havior of the DNWR algorithm with multiple subdomains: 


dttu — dxxW, 

u(x, 0) = 0, u t (x, 0) = 0, 
zi(0, t) = t 2 , u( 5, t) = t 2 e _t , 


x £ (0,5), t > 0, 
0 < x < 5, 
t > 0, 


(7.1) 
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• 3 subs error 
-•■3 subs bound 


.4 subs error 


—4 subs bound 

“-X ★ X \ 

★5 subs error 

\\ ★ \ , \ 

-5 subs bound 

★ \ _ \ 

• 6 subs error 


—6 subs bound 


v \ 



Fig. 7.1: Convergence estimates of DNWR for 9 = 1/2, on the left for various values 
of T for five subdomains, and on the right for various number of subdomains for T = 2 




Fig. 7.2: Arrangement A3: convergence of DNWR with various values of 9 for T = 5 
on the left, and for various lengths T of the time window and 9 = 1/2 on the right 


which is discretized using centered finite differences in both space and time on a 
grid with Aa; = At = 2xl0 -2 . We take the initial guesses g°j{t) = t 2 ,t £ (0,T] 
for 1 < j < 4, and consider a decomposition of (0,5) into five unequal subdomains, 
whose widths hi are 1, 0.5,1.5,1,1 respectively, so that h m i n = 0.5. On the left panel 
of Figure 7.2, we show the convergence for different values of the parameter 9 for 
T = 5, and on the right panel the error curves for the best parameter 9 = 1/2 for 
different time window length T. These convergence curves justify our convergence 
result for the arrangement A3 in Theorem 3.1. We observe two-step convergence for 
9 = 1/2 for a sufficiently small time window T. Coincidentally we observe exactly 
the same convergence behavior for other two corresponding arrangements Al and A2 
from Subsection 2.1, see Figure 7.3 and Figure 7.4 respectively. 

Next we show an experiment for the DNWR algorithm in two dimension for the 
following model problem 


d tt u— {d xx u + dyyU ) = 0,u(x,y,0) = xy(x- l)(y - n)(5x~ 2){Ax - 3), u t {x, y, 0) = 0, 

with homogeneous Dirichlet boundary conditions. We discretize the wave equation 
using the centered finite difference in both space and time (Leapfrog scheme) on a 
grid with Ax = 5xl0 -2 , Ay = 16xl0 _2 ,At = 4xl0 -2 . We decompose our domain 
D := (0,1) x (0,7r) into three non-overlapping subdomains fR = (0, 2/5) x (0, n), D 2 = 
(2/5, 3/4) x (0, 7r), fl 3 = (3/4, 1) x (0, n). As initial guesses, we take w®(y,t) = tsm(y). 
In Figure 7.5 we plot the convergence curves for different values of the parameter 9 for 
T = 2 on the left panel, and on the right the results for the best parameter 9 = 1/2 
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Fig. 7.3: Arrangement Al: convergence of DNWR with various values of 9 for T = 5 
on the left, and for various lengths T of the time window and 9 — 1/2 on the right 




Fig. 7.4: Arrangement A2: convergence of DNWR with various values of 9 for T = 5 
on the left, and for various lengths T of the time window and 9 = 1/2 on the right 


for different time window length T. 

We now compare in Figure 7.6 the performance of the DNWR algorithm with 
its counterpart NNWR method [ 4 , ] and the SWR algorithms with and without 
overlap [ 7 , ]. Here we consider the model problem 


d tt u - {d xx u + dyyu) = 0 ,u(x,y, 0) = 0 = u t (x,y , 0), 


with Dirichlet boundary conditions u(0,y,t) = t 2 sin(y), it(l, y 1 1) = y(y — 7r)t 3 and 
u(x,0,t) = 0 = u(x, 7r,f). We decompose our domain Q := (0,1) x (0,7r) for the two 
subdomains experiment into fli = (0, 3/5) x (0, n) and f h = (3/5,1) x (0,7r), and for 
the three subdomains experiment into fli = (0,2/5) x (0,7r), f l 2 = (2/5, 3/4) x (0,7r), 
fl 3 = (3/4, 1) x (0,7r). We take a random initial guess to start the iteration, and for 
the overlapping SWR we use an overlap of length 2Aa; in all the experiments. We 
implement first order methods with one parameter in optimized SWR iterations; for 
more details see [ ]. On the left panel of Figure 7.6 we plot the comparison curves for 
two subdomains, and the same for three subdomains on the right. 

Here we consider a comparison of performances between the DNWR and the 
NNWR algorithms for the wave equation. Table 7.1 gives a summary of the theoretical 
results from Section 3 and 6 and [ 1, 9], to indicate the maximum number of iterations 
needed for the ID and 2D wave equation to converge to the exact solution. 

Next we show a numerical experiment for the DNWR algorithm with different 
time grids for different subdomains and discontinuous wave speed across interfaces. 
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Fig. 7.5: Convergence of DNWR in 2D: curves for different values of 9 for T = 2 on 
the left, and for various time lengths T and 9 = 1/2 on the right 



iteration iteration 


Fig. 7.6: Comparison of DNWR, NNWR, and SWR for T = 2 in 2D for two subdo¬ 
mains on the left, and three subdomains on the right 


We consider the model problem 


d tt u - c 2 d xx u = 0, u(x, 0) = 0 = u t (x, 0), 

with Dirichlet boundary conditions u(0,t) = t 2 ,u(6,t) = t 3 . Suppose the spatial 
domain D := (0,6) is decomposed into three equal subdomains f \,i = 1,2,3, and 
the random initial guesses are used to start the DNWR iteration. For the spatial 
discretization, we take a uniform mesh with size Aa; = 1x10 , and for the time 

discretization, we use non-uniform time grids Ati,i = 1,2,3, as given in Table 7.2. 
For the non-uniform mesh grid, boundary data is transmitted from one subdomain to 
a neighboring subdomain by introducing a suitable time projection. For two dimen¬ 
sional problems, the interface is one dimensional. Using ideas of merge sort one can 
compute the projection with linear cost, see [ ] and the references therein. In Figure 
7.7(a) we show the non-uniform time steps for different subdomains. Figure 7.7(b), 
(c) and (d) give the three-step convergence of the DNWR algorithm for T = 2. 

8. Conclusions . We defined the DNWR algorithm for multiple subdomains for 
parabolic and hyperbolic problems, and analyzed its convergence properties for one 
dimensional heat and wave equations. We proved using numerical experiments that 
for a particular choice of the relaxation parameter, 9 = 1/2, superlinear convergence 
can be obtained for heat equation, whereas we showed finite step convergence for 
wave equation. In fact, the algorithm can be used as a two-step method for the wave 
equation, choosing the time window lengh T small enough. We have also extended the 
DNWR algorithm for 2D wave equation, and analyzed its convergence properties. We 
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Table 7.1: Comparison of steps needed for convergence of DNWR and NNWR for the 
wave equation. 


Methods 

2 subdomains, ID 

Many subdomains, ID 

Many subdomains, 2D 

DNWR 

T < 2khm\ n /c 

T < kh min /c 

T < kh m i n /c 

NNWR 

T < 4fc/i min /c 

T < 2kh min /c 

T < 2 kh min /c 


Table 7.2: Propagation speed and time steps for different subdomains. 




D2 

D3 

wave speed c 

1/4 

2 

1/2 

time grids At* 


39 x icr a 

1 x 10" 1 



(c) DNWR: 2nd iteration (d) DNWR: 3rd iteration 


Fig. 7.7: Convergence of the DNWR method applied to the wave equation for non- 
uniform time steps for 9 = 1/2 for T = 2 
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have also shown using numerical experiments that among DNWR and NNWR, the 
second converges faster. But in comparison to DNWR, the NNWR has to solve twice 
the number of subproblems (once for Dirichlet subproblems, and once for Neumann 
subproblems) on each subdomain at each iteration. Therefore the computational cost 
is almost double for the NNWR than for the DNWR algorithm at each step. However, 
we get better convergence behavior with the NNWR in terms of iteration numbers. 
Finally we presented a comparison of performences between the DNWR, NNWR and 
Schwarz WR methods, and showed that the DNWR and NNWR converge faster than 
optimized SWR at least for higher dimensions. 
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